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Abstract. Inspired by the penalization of the domain approach of Lions & 
Sznitman, we give a sense to Neumann and obUque derivatives boundary value 
problems for nonlocal, possibly degenerate elliptic equations. Two different 
cases are considered: (i) homogeneous Neumann boundary conditions in con- 
vex, possibly non-smooth and unbounded domains, and (ii) general oblique 
derivatives boundary conditions in smooth, bounded, and possibly non-convex 
domains. In each case we give apropriate definitions of viscosity solutions and 
prove uniqueness of solutions of the corresponding boundary value problems. 
We prove that these boundary value problems arise in the penalization of the 
domain limit from whole space problems and obtain as a corollary the existence 
of solutions of these problems. 



1. Introduction 

Inspired by the penalization of the domain approach of Lions & Sznitman [20] 
(see also [5T1 [53]), we give a sense to Neumann and oblique derivatives boundary 
value problems for nonlocal degenerate elliptic partial integro-differential equations 
(PIDEs in short). Because of the nonlocal nature of our PIDEs posed in a domain 
il, the boundary conditions have then to be imposed not only at the boundary 9f2, 
but possibly in all of the complement At least boundary conditions must be 
imposed in the union of the supports of the jump measures (see below). 

To be more specific, we consider PIDEs of the form 

F{x,u,Du,D'^u,I[u]{x)) ^ in n, (1.1) 

with extended Neumann/oblique derivatives boundary conditions 

Du{x)--/ix) = g{x) in . (1.2) 

Here $7 is a domain in and is a real-valued, continuous function defined on 
X K X X §^ X K, where is the space of TV x symmetric matrices. We 
assume that F is degenerate elliptic which, in this nonlocal setting, means that for 
any x e M^, u e M, p e M^, M,N eS^ and h,l2 G M, 

F{x,u,p,MJi) < F{x,u,p,N,l2) when M > TV, ?i > /a- 

Our assumptions cover the cases of general linear and nonlinear equations and, in 
particular, Bellman-Isaacs equations of control and game theory. 
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The equation is nonlocal because of its dependence on the nonlocal operator T[u\ 
which we assume to be of Levy-Ito type. For any smooth bounded function (f> and 
for any x d fl, 



I^x)^ [<|){x+J{x,z))~(|){x)^D(t>{x)■j{x,z)lBiz)]d^l{z), (1.3) 

where B C is the unit ball, 1b the indicator function of B, j is a Borel function, 
and fi - the Levy measure - is a positive Radon measure on M*^ \ {0}, and there 
exists c(j) > such that 

|zp A 1 < CX3 and \j{x, z)\ < c{j)\z\ for any |z| < 1, x e H. (1.4) 



A Taylor expansion shows that is well-defined under (1.4 1. We will assume 



without mention that (1.4) holds throughout this paper. The operator X is the 
generator of a stochastic jump process which solves a stochastic diffential equation 
involving a general jump term/Poisson random measure, cf. [141 124| . In particular, 
all generators of pure jump Levy processes [Tj are included. A typical example is 
the fractional Laplacian I = for a G (0,2), the generator of the symmetric 
a-stable processes, where 

dz 

j{x, z) = z for any x and z and fi{dz) - 



Most Levy models arising in Finance also fall into our framework, see e.g. |11) . 



For the boundary condition (1.2), we assume that 

(BCl) The functions 7 : ~> and g : ^ M are bounded Lipschitz 
continuous functions, and there exists > such that 7(0:) • n{x) > v for 
any x € d^. 

(BC2) For any x e fj'^, :— inft>o{Xj;(t) G 17} < +00, where Xx{-) solves 

X^{Q)^x and X^{t) = --i{X^{t)) for t>Q. (1.5) 



Assumption (BCl) is sufficient for (1.2) to really play the role of a boundary 



condition in the case of local equations. Assumption (BC2) states that integral 
curves of the vector field —7 starting from any point a; G il*^, will reach the bound- 
ary 9ri in finite time. This is a natural condition for a Neumann type boundary 
condition, and it is very closely related to the idea of the "penalization of the do- 
main" method of Lions & Sznitman [20] (see also [5T1[53]). This method is based 
on the observation that in the limit k 0, the vector field —^7 instantaneously 
returns the underlying stochastic process to Q, after an outside jump, and this is 
where (1.5) plays a role. We refer to Section [s] for more details in this direction. 



As in Lions & Snizman |20) . we use the notion of viscosity solutions. For nonlo- 
cal equations posed in full space, we refer to [6] (see also [3l[T9l[24]) and references 
therein for an account of this theory. A nonlocal Dirichlet problem was consid- 
ered in |S], where boundary conditions are given in all of VL'^ in an analogous way 
as in this paper. Here we consider two different cases: (i) a very general class 
of equations with homogeneous Neumann boundary conditions posed in convex, 
possibly non-smooth and unbounded domains, and (ii) a less general class of equa- 
tions with general oblique derivatives boundary conditions in smooth, bounded, 
and possibly non-convex domains. In each case we give appropriate definitions of 
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viscosity solutions and prove uniqueness theorems for the corresponding boundary 
value problems. Here we want to point out that the extended oblique derivative 
condition (1.2 1 influences the behavior of the solutions at infinity and therefore in- 
terferes in the conditions which are needed to have a well-defined nonlocal operator; 
we discuss this point at the end of Section [2j We also show that our formulation 
follows from a sequence of problems posed in the whole space obatined from the 
penalization of the domain method in a similar way as in [20j. As a consequence 
we also get some existence results for our problems. 

In some cases covered in this paper, a probabilistic description of the reflection 
problems based on stochastic differential equations can be found in [52 ■ Elsewhere 
in the litterature similar problems have been investigated for Levy operators where 
the measure ^.^ forces the underlaying process to stay in the domain either by a 
"smooth" restriction of its support or by "killing" all jumps leaving ^l. In these 
cases a Neumann boundary condition can be imposed only at the boundary just 
as for local problems. The first type of problems is considered in see also the 
book [13], and the killing approach is linked to the a-censored process [9] and the 
regional fractional Laplacian [151 [171 [16] . In these approaches, changes the domain 
lead to changes in the nonlocal operators and equations as well. E.g. equations with 
fractional laplacian operators can no longer be used directly. In our approach we 
can use the original (full space) operators at the cost of having boundary conditions 
in Qf^ . We also mention that in f4j, the authors along with E. Chasseigne investigate 
different ways of understanding Neumann boundary value problems for Levy-type 
nonlocal equations posed in the half space = M^^^ x M+. Four different ways 
of understanding the homogeneous Neumann condition are discussed, and the one 
involving normal projection of outside jumps is linked to this paper. 

In the case where the underlying process is a symmetric a-stable processes (a 
subordinated Brownian motion), the above mentioned approaches follows after a 
"reflection" on the boundary: The processes can be constructed from a Brownian 
motion by flrst subordinating it and then reflecting it. Another possible way to 
construct a "reflected" process it to first refiect the Brownian motion and then 
subordinate the refiected process. This approach is related to Dirichlet-Neumann 
operator, and it have been described e.g. by Hsu [18] using probabilistic methods 
and by Caffarelli and Silvestre by analytic PIDE methods. Especially the ideas 
of |10| have been used by many authors since. 

Our paper is organized as follows: Section [2] is devoted to a key technical lemma 
which allows us to control the solutions outside the domain Vt. We also discuss 
the connections between the extended oblique derivative condition (1.2), the be- 
havior of the solutions at infinity, and the conditions which are needed to have a 
well-defined nonlocal operator. In Section [3| we will focus on convex possibly non- 
smooth and unbounded domains but restrict ourselves to homogeneous boundary 
conditions. We define the concept of viscosity solution and proof a comparison 
theorem. The key argument here is to obtain by convexity a contraction property 
that force maximum points of the test function to be in $7. After this, the proof 
can be concluded in the standard (full space) way. In Section [4j we prove a com- 
parison theorem in the case of general oblique derivative conditions and smooth 
bounded possibly nonconvex domains. The proof uses the complicated test func- 
tion constructed by G. Barles in [5], along with the technical lemma of Section [2] 
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Finally, Section [5] is devoted to the analysis of an asymptotic result - the penal- 
ization of the domain method introduced by Lions and Sznitman. We prove that 
the above boundary value problems arise in the penalization of the domain limit of 
whole space problems and obtain as a corollary existence results for our Neumann 
problems. 

2. Preliminary Results 

In this section we state and prove two lemmas which play key roles in the proofs 
in the next sections. We recall from (BC2) that Ty := mitya{Xy{t) e H.} for y E . 

Lemma 2.1. Assume (BCl) and (BC2). 

(a) If u is a locally bounded, use funetion satisfying Du{x) ■ '){x) < g{x) in ft in 
the viscosity sense, then, for any y E fl and for any t < Ty, we have 

u{y)< f g{Xy{s))ds + u{Xy{t)). (2.1) 
Jo 

(h) If V is a locally bounded, Isc function satisfying Dv{x) ■ j{x) > g{x) in fj"^ in 
the viscosity sense, then, for any y E fl and for any t < Ty, we have 

v{y)> f g{Xy{s))ds + v{Xy{t)). (2.2) 
Jo 




Proof. The proof is inspired by an argument of G. Barles, S. Mirrahimi, B. Perthame 
and P.E. Souganidis 8.J. We only prove (a) since the proof of (b) is similar. Let 

t £ (0, Ty] and define 

w{s) := u{Xy{t - s)) for s > 0. 
If we can prove that the use function w is a subsolution of 
dw , 
ds 

then by the comparison principle 



-{s) ^ g{Xy{t^ s)) fors>0, (2.3) 



w{s) < w{0) + / g{Xy{t - T))dT for s > 0, 



since the right-hand side is the C^-solution of ( |2.3| ) with same initial data ^(0). 
Now (2.1) follows by choosing s = t. 
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To prove that w is a supersolution of (2.3), we take any smooth test-function 
and any point s > such that w — ip has a local strict maximum point at s. Note 
that Xy(t — s) € f2 since s > 0. Next we introduce the functions 



(x, s) = 



\Xy{t-S)-X\ 



+ <y9(s) and il)(x^s) — u{x) — 4>(x^s). 



For £ > small enough, classical arguments show that -0 has a maximum point 
near [Xyit — s),s) (depending on e) that we also call {x, s). Moreover 



\Xyit - S) - X\^^ 



— > and s — 



as 



£ 0. 



(2.4) 



Therefore x € Q for £ small enough, and since u is a subsolution of (1.2) and 
u — 0(-, t) has a local maximum at x. 



2{Xy{t-s)-x) 



lix) < g{x). 



Since r i— >■ ij^^x, t) is a C^-function having a local maximum at t = s, we also have 

2(Xy(t-s]-x) ■ , , , di}) , , „ 

^2 -Xyit-s)- Vi(=s) = ^(x, s) = 0, 

and we can conclude using Xy{t — s) = —j{Xy(t — s)) and Lipschitz continuity of 
7 that 



(p{s) 



< 



2{Xy{t ~ s)~x) 



2{Xy{t-s)-x) 2{Xy{t-s)-x) 



l{Xy{t-S)) 



{^{Xy{t-S))-^{X)) 



< 9{x) + 



\Xy{t ~ S) ~ X\' 



In view of (2.4), we can send £ — > to find that ip{s) < g{Xy{t — s)). 



(2.5) 



□ 



To allow for convex domains with corners and 7 = n, we need to relax the 
Lipschitz assumption on 7 in (BCl) and impose only a one-sided Lipschitz condition 



(BCl') The functions 7 : VL 
functions, 



and g : f2 



are bounded continuous 



- 7(y)) ■{x-y)> ~K\x ~ 2/1" 



(2.6) 



for some constant K and for all x, y G , and there exists v > Q such that 
7(x) • n[x) > V for any x G dO.. 



We state a slight generalization of Lemma 2.1 



Lemma 2.2. The results of Lemma 2.1 remain valid if we replace (BCl) by (BCl '). 

Proof. We first remark that (BCl') ensures the existence and uniqueness of the 
trajectory Xy as long as it remains in VL , existence follows from Peano's Theorem 
while (2.61 provides the uniqueness. Therefore the trajectory exists on [0,Tj^) and 
can be extended by continuity to t = r„. 
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To complete the proof, we let t < Ty and redo the proof of Lemma 2.1 The 



computations are exactly the same, e.g. to obtain (2.5), (BCl') is sufficient. To 
extend the result to t = Ty, we just remark that 



liiasup u{Xy{t)) < u{Xy{Ty)) , 

ttTy 



by the upper-semicontinuity of u. 



□ 



We conc lude this section by an i mpor tant discussion on the consequences of 



Lemma 



2.1 



If M is a solution of (1.1)-(1.2), then Du{x) •7(2;) — g{x) in and we 



have, for all y € 17 



u{y)= / g{Xy{s))ds + u{Xy{Ty)). 
Jo 



If (7 = on , then u can be bounded under suitable assumptions on the (other) 
data. This is the case we face in Section [3] below for convex domains and extended 
homogeneous Neumann boundary conditions. But if g is not identically on , 
then u can be unbounded and its growth is governed by the properties of g and Ty. 

The behavior of u at infinity is important in our framework to insure that the 
nonlocal terms are well-defined: We need some integrability property like e.g., for 
any a; G f2 and 6 > 0, 



\z\>S 



\u{x + i{x, z))\ d^{z) < +00. 



(2.7) 



Such condition now connects the assumptions we have to place on t^, g, j and /i. 



To fix ideas, we are going to assume in Section 5.2 that: 

(BC3) Either the function g has a compact support in or there exists c > 
such that for from (BC2), 



Tx < C(l + \x\ 



and 



sup 



\z\>S 



\j{x, z)\ d^j,{z) < 00 for any S > . (2.8) 



We briefly comment on this assumption. When g has compact support, the 



solutions are expected to be bounded by Lemma 2.1 and no additional assumption 
on j and fi is needed. On the contrary, if e.g. g = 1, then the integral of g in (2.1 ) 



suggests that u and Tx have the same growth and (BC3) imposes a linear growth. 
Next one has to impose suitable hypothesis on j and /i to satisfy (2.7 1. This is 



obtained through the second part of (BC3) on the /i-integrability of j away from 0. 
3. The homogeneous Neumann condition in convex non-smooth 

DOMAINS 

In this section we consider the homogeneous Neumann problem, namely equation 



( 1.1 ) and boundary condition ( 1.2 ) with g = and 7 = n, the unit outward normal 

(3.1) 



vectorfield in fl'^ (see below) 



Du{x) ■ n{x) = in Q", 

in the case when is a convex, possibly unbounded and non-smooth domain. 
At a; G dft, the set of outward normals Nq{x) can be defined as 

Nn{x) := |n e : |n| = l,n • (x - y) > for aU y e Hj. 
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This set is a singleton at any point where dfl is and part of a convex cone where 
dfl has a corner. Let d be the distance function to fl, and note that d = in f2 and 
d > in Q . Moreover, d is convex and belongs to C°(M^) n C^n'') since f2 is a 
closed convex subset of M^. In we now define the outward unit normal vector 
n in the only natural way by setting n — Dd. Note that the two definitions are 
consistent in the sense that 



'X^) = { 



n G 



lim Dd{xk) for some Xk x> for all x € Sfi, 



and since d is convex, the function n satisfies (BCl') with K = 0. 

To define the concept of viscosity solutions for this problem, we need the opera- 
tors Is, , and T defined as follows 



ismx) 

I'[u]{x) 



{x + j{x, z)) - (t>{x) - D(j>{x) ■ j{x, z)Ib{z) dfi{z), 



\z\<S 



u{x + j{x, z)) — u{x) — D(f>{x) ■ j{x, z)1b{z) dfi{z), 



\z\>S 



T[u,(l3]{x) = F{x,u{x),D(j){x),D^(j){x),l5[(l3]{x) +l'^[u]{x)). 



Under assumption ( 1.4 ), [0] is well-defined for g C^. For [u] we need some 
integrability condition a la (2.7 1, see the discussion at the end of Section [21 In this 

C II'' ' 

section 5 = on J7 , so (2.7) will be automatically satisfied whenever u is bounded. 



Definition 3.1. (i) A locally bounded, use function u : — !■ M is a viscosity 
subsolution of (1.1|-(1.2) if it satisfies (2.7), and for any test function (f) £ C^(M^) 
and for any maximum point G M.^ of u ~ (j) in Bc(j)s(xo) where c{j) is defined 
in ( 1.4), we have 



■r[u,^]{xo) < 



if xq e 



min T[u, (f)]{xo), inf D(f){xQ) • n < if xq £ dfl, 

\ nGNiiixo) J 

D(I){xq) ■ n{xQ) < if Xq G Q 

(ii) A locally hounded, Isc function v : M.^ — > M is a viscosity supersolution of 
(1.1)-(1.2) if it satisfies (2.7), and for any test function (j) € C^(M^) and for any 
minimum point xq G of the function u — (j) in B^(^j-^g{xQ) where c{j) is defined 
in ( 1.4), we have 



T[u,4>]{xa) > 



if Xq e 



max ( T[u, (f)]{xo), sup D(f){xQ) • n ) > if xq G dfl, 

V neNnixo) 

D4>[xq) ■ n(xa) > 



if Xq £ 



(Hi) A viscosity solution u of (1.1 1 -(1.2) is a locally bounded function whose upper 
and lower semicontinuous envelopes are respectively sub- and supersolution of the 
problem. 

This definition is a natural extension of the definition given in 5j to the Neumann 
type boundary value problem. 
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Remark 3.2. Two useful equivalent definitions can be given: (1) We can replace 
I^[u] by in the above definition if local maximum/minimum points are re- 

placed by global ones. (2) In the subsolution definition, {D(t){xo), D'^(j>{xo)) can be 
replaced by elements {p,X) in the so-called super-jet J^u{xo) if X < £'^0(a;o). 
In the definition of supersolutions, you can similarly use {q,Y) G J^u(xo) if 
Y > D'^(j){xo). The second definition is useful for comparison proofs, and the 
proofs of these claims easily follow from the arguments for similar results in [5] . 

We now state the assumptions - remarking that the assumptions on F will be 
as general as for the whole space case fl = M.^ without boundary conditions. For 
convenience we use the assumptions of [6 , but see Remark 3.6 below. For the 
nonlocal part we assume that 

(Al) There is a constant c > such that for all x,y ^ M^, 



\j{x,z)\'^n(dz) 



'j{y,z)\^ 



y\ 



\j{x,z) - j{v,z)\ 



'\B 



v\ 



li{dz) < c. 



The non-linearity F is continuous and satisfies the following classical assumptions 



(A2) There exists a > such that for any x € 



u,v eR, p eR'^, X and 



(A3-i; 



F{x,u,p, X,l) ~ F{x,v,p, X,l) > a{u ~ v) when u>v. 

For any i? > 0, there exist moduli of continuity uj^lor such that, for any 
|a;|, \y\ < R, \v\ < R, I (^M. and for any X, F G satisfying 



X 








Y 



< 



1 



(3.2) 

as /3 for 



(3.3) 



/ -/ 
-/ / 

for some e > and r(j3) a.s f3 0, then, if Si{/3) 
i = 1 and 2, we have 

f (y, V, s~'ix -y) + Slip), Y, I) ~ F(x, v, e~\x - y) + S2{l3),X, I) 

< uj{f3) + ujr{\x ^ y\ + e-^\x ~ yf). 

(A3-2) For any i? > 0, F is uniformly continuous on M" x [-R, R]x BrX Dr x M 
where Dr := {X E \X\ < R} and there exist a modulus of continuity 
ujR such that, for any x,y € M'', |u| < R,l £ R and for any X,Y e S'^ 
satisfying ( |3.2[ ) and e > 0, we have 

F{y,v,e-\x-y),Y,l)-F{x,v,e-\x-y),X,l) <ujRie-^\x-y\^ + \x-y\). (3.4) 

(A4) F{x,u,p, X,l) is nondecreasing and Lipschitz continuous in I, uniformly 
with respect to all the other variables. 

(A5) Mf :=sup,gf2|F(a;,0,0,0,0)| <oo. 

Assumption (A3-1) and (A3-2) are two versions of assumption (3.14) in the Users' 
Guide [12 for possibly unbounded domains. These assumptions along with (A4) 
imply that equation (1.1) is degenerate elliptic. Assumptions (A3-1) allows more 
general x-dependence in the equation (e.g. general x-depending HJB equations), 
while (A3-2) allows more general gradient dependence in the equation (e.g. non- 
x-depending equations with superlinear gradient terms). In the local case with no 
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I-dependence in the equation, assumptions (A2) - (A5) imply comparison, unique- 
ness, and existence (via Perron's method) of a bounded viscosity solution of 



( 1.2 ), cf. e.g. [T^. In the nonlocal case when f2 = M (and no Neumann conditions. 



n'^ = 0), we have the following rather classical result which we will need later. 

Proposition 3.3 (Results for n ^ M^). Assume 17 = and (Al), (A2), (A4) 
hold along with either (A3-1) or (A3-2). 

(a) If u and v are respectively an use hounded above subsolution and a Isc bounded 
below supersolution of (1.1 ) in ^ M.^ , then u < v in M.^ . 



(b) Assume also (A5) holds, then there exists a unique bounded viscosity solution 
u of (1.1) in ri = M.^ satisfying 

Hx)\ <— m M^. (3.5) 
a 

Part (a) was proved in [S] (see Section 5), and Part (b) follows from part (a) 



and Perron's method since Mp/a and —Mp/a are super and subsolutions of ( 1.1 1. 
Similar results have been given e.g. in [3} [24 l [T9 t l5] . 

Now we come to the first main result of this paper, a comparison result for the 



boundary value problem ( 1.1 )-(^3.1h. 



Theorem 3.4 (Comparison I). Assume (Al), (A2), (A4) hold along with either 
(A3-1) or (A3-2). If u and v are respectively an use bounded above subsolution and 
a Isc bounded below supersolution of ( |1.1[ )-( [33| ), then u < v in M^. 

Uniqueness of solutions follow, and since are sub/super solutions of (1.1) 

when (A5) holds, we also get L°°-bounds. 

Corollary 3.5. Assume (Al), (A2), (A4) hold along with either (A3-1) or (A3-2). 



(a) There is not more than one bounded solution of (1.1) 



(b) If also (A5) holds, then any solution u of (1.1) satisfies (3.5) 



Remark 3.6. Under assumption (A3-1), the above results also holds if assumption 
(Al) is replaced by the much more general assumption 

(Al-2) The function x i— > j{x, z) is continuous for a.e. z and there exists a constant 
c > such that 



\i{x,z) - j{y,z)\^^i{dz) < c\x - yp. 

IB 

The proof in Section 5 in |S] can be modified easily to cover this case by a clever 
trick which can be found e.g. in Section 6 in |19| . As opposed to assumption (Al), 



assumption (Al-2) covers the fractional Laplace case. If we also relax (1.4) so that 
the constant c{j) is finite only for compact subsets oi x G ft, then the above results 
also cover the case when j has linear growth in x. 

Proof of Theorem \3.4\ Wc introduce the following initial value problem (cf. (BC2)), 
Xyit) ^ -n{Xy{t)) for t>0, Xy{Q) = y. (3.6) 
Note that the projection on the closed convex set fi, P : fi^ — >■ dVl, is also given by 

Py = Xy{Ty) ATiy ?/ £ 17 , 
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Figure 2. The curves of (BCl) with ^ ~ n starting at yi and y2 
and ending at a corner point of 9f2. 



where we recall that 

Ty = inf{t > : Xy(t) e dn}. 

Since is convex and \n\ = 1, it follows that {Xy{-)}y defines a family of constant 
speed, finite length, and non-intersecting paths in o!^ having the form 

Xy{t) = y - tn{y) for t(^[Q,Ty\. (3.7) 

Obviously Ty < oo for all y so that (BC2) is trivially satisfied when 7 = n. 
We argue by contradiction assuming that 

M :— sup{M(a;) — v{x)} > 0. 



Since n satisfies (BCl'), Lemma 2.2 applies with g = and we find that u{x) — 

v{x) < u{Px) — v{Px) in ^ , and hence that M = niax{u(a::) — v{x)}. 

a 

Let X : M ^ M be a bounded smooth function such that x(^) = for ^ < 0, 
x'{t) > for t > 0, and x{t) > 2(||m||oo + ||w||oo) for t>l. We double the variables, 
introducing the function 

'fix, y) - u{x) - v{y) - tl^! _ _ ^^.p + 1)) _ _ ^ 1))^ 

where £,/3 > 0, and xq is any given point in fJ. It is easy to see that, for /3 
small enough, M^^^ = maxK2N "f{x,y) exists and is attained at some point {x,y) e 
X (that depends on e and (3). The crucial and new step in the proof is 
to show that (a;, y) e f2 x f2. If this was not the case, then two applications of 
Lemma [2^ yields that 



u(x)-v{y)<u{P^)-v{Py). (3.8) 
Moreover, since £7 is convex and a;o G il, 

\x - y\>\Px - Py\, \x - xq\>\Ps - xq\, and \y - xo\ > \Py - xo\, (3.9) 
and then, for /? small enough, we have the contradiction 

M,,^ = *(x,y)<*(Ps,Pj^). (3.10) 
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Since x,y & il, the rest of the proof fohows classical arguments. Assume x G 
and let 

0(x, y) = + x(/3(|x - xop + 1)) + - xoP + 1)). 

Note that by convexity of ft, 

{x- y) ■n>0 for all y e ^, n e Nn{x). 

Moreover, this inequality is strict ii y ^ il. Finally, since x'(i) > for i > 0, we 
use the fact that xq E fl to find that 
2 

DxHx^y) ■ n = -^{x - y) • 7i + x'(/3(|S - a;oP + l))2/3(x - xo) ■ n 
> for all n G Nq{x) and /3 > 0. 



Therefore, from Definition 3.1 the equation has to hold at i, i.e. J-[u, y)]{x) < 
0. A similar argument shows that J^[v, —(j){x, ■)]{y) > if y € 

Now we are in the situation that x,y € ^ and that the equation is satisfied at 
these points. The conclusion of the proof is then exactly as for the M.^ case, and 
we omit the standard details. Under the present assumptions, essentially all the 
remaining details can be found in Section 5 in 1^ . But see also [31 [121 [Mj for very 
similar results. □ 

Remark 3.7. The key ingredients of the above proof are 



(i) Inequality (3.8) that comes from Lemma 2.1 or 2.2 and that allow us to 



compare values of u and v outside ft with those on inside 17. 



(ii) Inequality (3.10) that comes from convexity and contraction properties (see 
(3.9 1). In the above proof, the contraction property of the projection on 
the closed, convex set f2 was playing the key role (allowing us to use a very 
simple test function), but in general the contraction property comes from 
the control on the Xy trajectories w.r.t. y. 

(iii) As in the classical Neumann/oblique derivatives boundary conditions cases, 
the test-function has to be build in order to allow us to "avoid" the bound- 
ary condition (cf. ( 3.11| )). 



These three ingredients are the same in any proof but with different arguments to 
handle them. We are going to focus on these arguments. 

Remark 3.8. If is bounded, we can relax assumptions (A3-1) and (A3-2) in the 
standard way and the comparison result will still hold. E.g. since we no longer 
need to prevent maximum points from escaping to infinity, we can set all functions 
r,si,S2 and uj equal zero in (A3-1). 

4. General oblique derivative conditions in non-convex smooth 

DOMAINS 

In this section we consider the general oblique derivative problem of the form 



( 1.1 )-( 1.2 ) on a bounded, possibly non-convex , C^-domain fi. Compared to section 
[3| the domain and boundary condition and are more general, but the class of 
equations (see below) and the boundary regularity are more restricted. 



of viscosity solutions 



Assuming that (1.4) and (BCl) hold, and we now have the following definition 
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Definition 4.1. (i) A locally bounded, use function u : 



is a viscosity 



subsolution of (1.1|-(1.2) if it satisfies (2.7), and for any test function cj) ^ C^(M^) 
and for any maximum point xq G M.^ of u — (j) in Bf.(^j^g(xQ) where c{j) is defined 
in Ol), 



J'[u,(f>]ixo) < 

mm I J^[u,(j)]{xQ),D<j){xo) ■ j{xo) 



if Xq e V,, 

j{xo) ) < if xoe an, 

if Xq G 



pN 



^D(I){xq) ■ 7(xo) < g{xQ) 
(ii) A locally hounded, Isc function v 

minimum point xq E M.^ of the function u — <j) in B^Q'^g{xQ) where c{j) is defined 
in Ol, 



R is a viscosity supersolution of 
(1.1)-(1.2) if it satisfies (2.7), and for any test function <j) £ C^(M^) and for any 



T[u,(f>]{xo) > if Xq en, 

ma.x (^J^[u, (j)]{xQ), D(I){xq) ■ j{xo) - g{xo)^ > if xq e dfl, 
Dcj){xQ) ■ 7(xo) > g{xQ) if xq G n" 



(Hi) A viscosity solution u of (1.1 1 -(1.2) is a locally hounded function whose upper 
and lower semicontinuous envelopes are respectively suh- and supersolution of the 
problem. 

To handle non-convex domains and more general boundary conditions, we will 
use a rather complicated test-function which is no longer only a function oix—y plus 
small terms. For the proofs to work out we therefore need to replace assumption 
(A3-1) and (A3-2) by a more restrictive assumption similar to the one used in the 
local case [2] 

(A3-3) For any R, K > 0, there exist moduli of continuity rn^i^K such that, for any 
x,y & il, \u\ < R, p,q £ M^, / G K, and matrices X,Y G satisfying 



y\ < 



K 



-Id< 



we have that 



X 






-Y 



\p-q\<K7jeil + \p\/\\q\), and 
K 



< 



Id -Id 
-Id Id 



Kr]Id, 



F(y, u, q, Y, I) - F{x, u,p, X, I) < tur^k [v + - 2/1(1 + \p\ V \q\) + '—^ 

We have the following comparison result. 

Theorem 4.2 (Comparison II). Assume (Al), (A2), (A3-3), (A4), (BCl), and 
(BC2) hold. If u and v are respectively a locally bounded use subsolution and a 
locally bounded Isc supersolution of ( |1.1[ ) -( 1.2 ), then u < v in . 



pN 



This result will be proved in the subsections below. We start by introducing the 
test function we need for the proof. 
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4.1. The test-function. As for local oblique derivative boundary conditions (see 
e.g. [2] and references therein), the proof of our comparison result requires a rather 
complicated test-function. Fortunately there are no major differences between the 
test-function for the local and nonlocal cases, and we now recall a few facts about 
the test-function of 2 and describe the adaptations we need to make here. 

We start by changing our definition of the "distance to the boundary" d. Now d 
will be a bounded function which is equal to the signed distance function to dfl in 
a neighborhood of dH. {d > mfl and d < inil ) and where n{x) := —Dd{x) ^ 
in Qf^. Note that n(x) is the outward unit normal vector to dVL for any x G d^. 
The test-function ^/^g^^ e C'^iM?'^) of ^2] can then be defined as follows, 



- e-^^W-)+'^('^)l fc^) {d{x) - d{v)) (4.1) 

^ ^-K, A{d{x)-- d{y)f _ ^^^^ ^^^^^ ^ ^^^^^ 



for parameters ?7, e > (small), constants A, Ki, K2 (large), and where the function 
Cjjs (see [2] page 214) is a suitable smooth approximation of a bounded Lipschitz 
extension of the solution t ~ C{x,p) of the equation 

7(0;) • {p + tn{x)) — g{x) = for p e M^, x near dVt. 

The key properties of the test-function are given in the Lemma below. 

Lemma 4.3. Assume (BCl) and let R> 0. //ry, e > are small enough, then for 
A,Ki,K2 large enough, then the function 4'e,ri defined in (4.1 1 has the following 
properties 

(i) For any x,y M.^ , 

^eAx,y) > ''"^ ~ Ks^ - K2f]e [d{x) + d{y)]. (4.2) 

(ii) For e,rj £ (0, 1) and \x — y\ < rye, 

iD^i^eA^^ y)\ + IDyi^eA^, y)\>~K + i^-^^-^, 

\Dx1pe,vi^^y)\ + l^yV'e,r)(a;,y)| < 

C^^+c(l + 772ifie2^ill'^ll- +e7yir2), (4.3) 
\D^^,^^ix, y) + Dyij,^^{x, y)\ < /iT *^"/* + K{rje + e^), and 



f^Id < 7^V.,,,(x, y)<^( ^'1^ // ) + Krjid. (4.4) 



(Hi) There is S > such that for \x ^ y\ < S and x,y in a neighborhood of dVl, 

7(x) • D^i>,^^{x,y)) > g{x) if d{x) < d{y) , (4.5) 
-7(2/) • Dy^,^r,{x,y)) < g{y) if d{y) < d{x) , (4.6) 
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and if in addition \x — y\ < r/e, then 

- ■ {D^tps^jj{x, y) + Dyips^^{x, y)) 



V 



(4.7) 



Except for (4.7), these estimates have essentially been proved in Section 5 in [2]. 
Some new features that only marginally changes the proofs are: (i) x,y can now 
belong to 51^, (ii) inequality (4.2) is slightly more accurate, and (iii) inequalities 
(4.5) and ( |4.6[ ) are now given in a neighborhood and not only at dQ. Moreover, the 
constants K will in general depend on Ki and K2, and the precise dependence is 
not important except for the term (4.3 ). The importance of this dependence is both 
new and central to this paper (cf. the proof of Lemma 4.4 a)). We will therefore 
prove both (4.3) and (4.7) here. 



Proof of (4.3) and (4.7). To simplify the computations, we write V'e,?; in the fol- 
lowing way 



V'e,r,(a;, y) = xix - y, d{x) - d{y), 



x + y 



dix)+diy)) 



where 

X{X,Y,Z,T) :=e 
In this notation 



v2 

£2 



Z, e 



-KiT 



2X 



Y 



1 



Dxipe,ri{x, y)^xx- xmix) + -xz - xmix), 

Dxipe.,r,{x,y) + Dyip^^r,{x,y) = ~XY{n{x) ~ nijj)) + Xz - XT{n{x) + n{y)). 

By the assumptions on 7 and g and the construction of Cj,e = Crieix,p) in [2 , there 
is a C > such that 



I a 



\DxC,r,\<C{l + \p\) and 



\DpC,^\ < G. 



Hence there are constants Ci and C2 such that 

'21X1 



\xx\ + \xy\ < Ci +6*26 
\xz\ < (C. + C^e 



-KiT 



21X1 



2(1 + A)|y| y 



1^1, 



IxtI <ifie-^^^C2( 



(i + A)y2 



K2T]e. 



Since \X\, \Y\ < C\x — y\, estimate (4.3) now follows. 

To prove (4.7), we note that by using Cauchy-Schwarz inequality on the DpCe-q- 
term and taking A large enough, 

'X^ 

t2 



XT 



-Kie 



-KiT 



< 



X2 



DpCne ■ ^2 



2X AY^ 

Y+ 



K2V£ 



AY' 



K2f]e. 



2 \ 

Let W = {x : dist(a::, dVl) < r}, and let r > be so small that 7 • n > | in W. Such 
a set exists by (BCl) and continuity of 7 and n. After an easy computation based 
on the above estimates, the Lipschitz continuity of n {\n{x) — n{y)\ ^ \X\), the 
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inequality 7 • n > |, Cauchy-Schwarz inequality, and finally, taking Ki,K2 large 
enough so that the XT-term dominates, we conclude that (4.7 1 holds in W. □ 



The next lemma plays a key role in the comparison proof. 



Lemma 4.4. Assume (BCl) and (BC2), let be defined in Lemma 2.1. and 
T \-nin{Tx, Ty). 

(a) For any K > 0, there are constants Ki,K2 large enough, such that for any 
£,77 > small enough, ifx,y € Q'^ are close enough to dVl and \x^y\ < rje/2, then 



(4.8) 



(b) For any rj > 0, there are constants Ki, K2 large enough, such that for any e > 0, 
if x,y € fi*^ are close enough to dfl and Ty < t^, then 



^ke,r,{X^{T.).Xy{T)) < i^,^^{X^{T),Xy{T)) 



g(X,{t))dt. 



(4.9) 



(c) For any rj > 0, there are constants Ki, K2 large enough, such that for any e > 0, 
if X £ Q'^ and y & ^ are close enough to dfl, then 

■4'e,n{^x{Tx),y) < i^e.7i{x,y) - / g{X^{t))dt. 

Jo 




Figure 3. Curves of (BCl) with different starting points in the 
oblique case. 



Proof Consider a neighborhood of dil, Wr = {x : dist{x,dfl) < r}, and let r > 
be so small that (4.71 holds, d{x) — ±dist{x,dil), and 7 • n > | in Wr- Such a set 



exists by the definition of d, (BCl), and continuity of 7 and n. In the set Wr H fl'^ 
the distance to boundary f{t) = dist{Xx{t),dQ) = —d{Xx{t)) is decreasing, 

V 



fit) = -Dd{X,{t)) ■ X,{t) = 7i{XS)) ■ (-7(X,(t))) < 



(4.10) 



and hence Xx{t) e n for aU t e [0, t^] and x eWrDil". 

Next we note that if L is the Lipschitz constant of 7, then by Gronwall's inequal- 
ity, 

\X,,{t)~Xy{t)\<e'^'\x-y\. (4.11) 



We estimate t^, and hence also Ty and r, by integrating (4.10) from t to Tx and 
noting that /(tj.) ~ 

-{tx -t) < f{t) = d\st{X^{t),dn) < dist(x,9f7) for t e [0,Tx]. 



16 



G. BARLES, C. A. GEORGELIN & E. R. JAKOBSEN 



Hence if r is small, t will also be small in Wr n fi^. In the rest of the proof we take 
x,y ^ Wr n and then we take r so small that also — Xy{t)\ < rje for all 

t e [0, r] and x,y eWrnfl" such that \x - y\ < ^. 

We now prove part (a). We start by using the definition of Xx{t) (see (BC2)) to 
show that 



dt 



We may use (4.3) (check!) and the Lipschitz continuity of 7 to have 



\Dy^P,,^ ■ [j{Xy{t)) ~ j{X^m\ 

■\X,it)-Xy{t)\ 



<L\X,{t)-Xyit)\-C(^^ 



< LC 



\X,it) - Xy{t)\^^ 



and by (4.7) we immediatly find that 

< _ ^ -Kild(X^{t))+d{Xym \^x{t) - Xy{t)\^ _ 

Since 7 • n > ^ , we then find that 



[i^,,,ix^{t),Xym 

_ '^-K,\d(Xjt))+d(X,(m\ \^x{t) - Xy{t)\ 



dt 

< [LC -Ki'-e 

+ {lc{i + er/X2 + Ty'ifie^^ill^ll-) - if2|)r?e 

< —Krjs 

for any given constant K since we can take first e, i] small enough and then Ki and 
finally K2 as large as we want. The conclusion follows by integrating from to r. 

To prove (b), we notice that t = Ty < t^. Since X = —'-f{X) and d{Xx{t)) < 
d{Xy{T)) ~ Q for T ~ Ty < t < Tx, we can use (4.5) to find that 

^ [AAXAt),Xy{T))] - -D^i^.^^-jiX^it)) < -g{Xx{t)). 
Part (b) now follows by integrating from Ty to r^.. The proof of (c) is just like the 

4.2. Proof of Theorem 



proof of (b) replacing Xy^r) by y and setting t 0. □ 

In order to show that u{x) — v{x) < in M^, we 



4.2 



first notice that, by (BC2) and Lemma 2.1 



u{x) - v{x) < u{Xx{t,j;)) - v{Xx{tx)) for any x efl , 

and hence since X^ir^) E dH,, it follows that u — u is bounded from above in 
and 

M — sup{u(x) — v{x)} = max{u(a;) — v{x)}. 
In the rest of the proof we argue by contradiction assuming that 



M > 0. 
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Then we define 

wpix) = u{x) - v{x) ~ 2x{-f3d{x)), 



where /? > (small), x is the function we introduced in the proof of Theorem 3.4 



and d is the signed distance function to {d < Oinfl ). Since the x-term vanishes 
on Q and is strictly positive on f2 , has maximum points only on Q and these 
points arc also maximum points of m — u. 

Now we double the variables introducing the function 

y) = u{x) - v{y) - il;^^^{x, y) - x{-l3d{x)) - x{~Pd{y))- 

By standard arguments involving the definition of x a-nd the properties of V'e,?; given 
in Lemma |4.3| (in particular (4.2 1), this function achieves its maximum at a point 
{x^y) e X (depending on e, rj and Moreover, for fixed rj and /3, 

^0 as e 0, 



and x^y converges (along subsequences) to a maximum point x of wp{x), i.e. to a 
point in H.. In particular, x, y will be arbitrarily close to dfl if e close enough to 0. 

We will show that x, y are in D, when e > is small enough. Again we argue 
by contradiction assuming that x, y are not both in f2. Assume e.g. that x,y G fl'^ 
and that Ty < t^. We will get a contradiction to the maximum point property by 
showing that 

^X,y)<^X^{T,),Xy(Ty)). 

To do this, we start by using Lemma [23] for both u and v to see that 



$(x,y) < uiX.^iT^)) ~ v{Xy{Ty)) 



+ 



" (g(Xs(5))-.g(X^(s)))ds+ / %(Xs(<))dt 

J Ty 

'0e,r,(S, 27) - x(-^d{x)) - x{-Pd{y)). 



But from Lemma 4.4 using first part (b) and then part (a), 

g{Xs,{t))dt < ^]^,r,{Xs;{Ty),Xy{Ty)) - tl)e,r,{Xs{Ts), Xy{Ty)) 

< ^e,r,{x,y) - 2KTT]e - ■4!^^jj{Xs.{Ts;),Xy{Ty)), 



and by Lipschitz regularity of g and 7 and the estimate (4.11 ), 

' [g{X,{s)) - g{Xyis)))ds < TyLgc'^-^^lx - y\. 

Hence we find that 
<^{x,y) 

< ^X^{T^),Xy{Ty)) - 2kTyf^e + TyLgC'^-^^ | ^ " ^ | " X ( " ^ ) ) " Xi-pCliy)), 

and since |a; — j/| < rye, we get the contradiction by choosing K large enough. A 
similar argument covers the case when Ty > r^, and we can conclude that at least 
one of X and y belongs to O. 
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Next we show that it is not possible that e.g. x & fl while y £ This time we 
use Lemma 12. 11 for only u to see that 







^x, y) < u{Xs{ts)) - v{y) + / g{X:,{t))dt - y) - x{-Pd{x)). 

Jo 

But by Lemma [4!4] |^c), 

9{Xx{t))dt - ips^jj{x,y) < -^pe,r,iXs(Ts),y), 

and hence we find again a contradiction 

<^{x,y)<'l>{X^{T^),y)-x{-Pdix)) <HX^{rs),y). 

The case that y G ^ while x G gives a contraction in a similar way, and in view 
of previous arguments we can conclude that x,y € fi, at least when e > is small 
enough. 



Since ipe,ri satisfies by (4.51 and (4.6 1, it follows that the equation (the sub and 
supersolution inequalites) , and not the boundary condition, has to hold if ai or y 
belongs to dil and hence for all x,y G Q. By assumption, u, v are bounded on f2 
so that assuption (A3-3) can be applied with R = max-^{\u\ + |ti|). At this point 
we can conclude the proof as in the M^-case, sending first e — >■ 0, then 77 — 0, and 
finally /3 — 0. We omit the standard details only noting that under the present 
assumptions, essentially all the remaining details can be found in Section 5 in [6]. 
But see also [El |3l [24] for very similar results. 

5. Penalization of the domain 

In this section we show that our way of defining Neumann type boundary condi- 
tions is consistent with the so-called penalization of the domain method introduced 
by Lions and Sznitman in [20] • We extend the results of [50] to our non-local set- 
ting, proving the convergence of a sequence of solutions of penalized M^-problems 
to the solution of (1.1). We give separate results in the convex case of Section [s] 



and the oblique case of Section [ij 

5.1. Neumann conditions on convex domains. In this section we assume that 
51 is convex and possibly unbounded. Let d be the distance to f2 defined in Sectionjs] 
and n = Dd in . Note that d = in and define d = mm{d, 1). By the Lipschitz 
continuity of d and the convexity of d, the continuous vector field x 1—^ d{x)n{x) 
(extended by to fi) satisfies ( [2.6[ ) in E^. This property will play a key role below. 

Moreover, we assume that (A1)-(A5) hold, and if necessary, we extend the data 
and F to M.^ in a way that preseverves these properties. We study the following 
equation for the penalization of the domain, cf. |20j 

F{x,u, Du, D'^u,I[u]{x)) H — d{x)n{x) ■ Du^O in M^. (5.1) 

K 



where < k <C 1. Since d{x)n{x) satisfies (2.61, Equation (5.1) with k > fixed 
satisfies (A1)-(A5) as long as F does. 

Theorem 5.1. Assume that (Al), (A2), (A4), (A5) hold along with either (A3- 



1) or (A3-2). Then the viscosity solution of (5.1) converge locally uniformly 
to a hounded continuous function u which is the viscosity solution of (jOJ -([L2]) 
according to Definition\3.1\ 
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Remark 5.2. This result provides an existence result for (1.1|-(1.2|. In contrast to 
the more difficult Dirichlet case in |S] , we have existence also when the there is loss 
of boundary conditions. 

We need the following auxilliary result that follows from Proposition |3.3[ 

Lemma 5.3. Assume that (Al), (A2), (A4), (A5) hold along with either (A3-1) or 
(A3-2). Then there exists a unique hounded viscosity solution of (5.1 1 satisfying 

\un{x)\ <— m M^. 
7 



Proof. Note that Mk is bounded uniformly in k, and that we may rewrite (5.1) in 
the following equivalent way 



Gi^{x,u,Du^D^u^I[u]{x)) =0 in 



i,N 



where 



G^{x,r,p,X,l) = 



F{x,r,p,X,l), for X £ Q, 

\^dl^Fi^'r,p,X,l) + n{x)p, for X en". 

Now we introduce the half relaxed limits 



(5.2) 



(5.3) 



/(a;) := liminf*/K(x) = liminf f^{y), f{x) := linisup*/K(2^) = limsup f^{y). 

y ^ ^ y — ^ X 

« K -J- 

Note that F = F and 

{F{x,r,p, X,l) when x € fl, 

min{F(x, r,p, X, /), MneN^tix) n ■ p} when x £ dQ, 

n{x)p when x £ , 

and in a similar way we find that G is like G with max/sup replacing the min/inf . 
As a consequence of the stability of viscosity solutions, see e.g. Theorem 1 in 6], 
u = lim sup* is a viscosity subsolution of 



G{x,u^Du,D^u^I[u]) = {i in 
while u = lim sup^ is a viscosity supersolution of 
G{x,u,Du,D^u^I[u]) = Q in 



By Definition 3.1 this means that u and u are sub- and supersolutions of ( 1.1 1-( 1.2 1, 



and hence by comparison, Theorem |3.4[ 

u < u. 

The opposite inequality is true by definition of u, and hence we have u = u =: u. 
It follows that u is continuous and — >■ u locally uniformly, as is standard in 
viscosity solution theory. □ 
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5.2. Oblique boundary value problems in bounded smooth domains. In 

this section, we assume as in Section |4j that i7 is a bounded domain. We study 
the fohowing equation for the penahzation of the domain, cf. [5D]: 



,T[uk\(x)) 



K 



Du^ - .g] = 



(5.4) 



where < k ^ 1 and d is defined as in the previous section. 

We want to prove that we can obtain the oblique boundary value problem ( 1.1 ) 
from the penalized problem (5.4) in the limit as /c — >■ 0. In (1.1) (Definition 4.1), 
only i^'s values at 51 play any role, and we may modify equation (5.4) in il"^ and 



still obtain (1.1) from (|5.4[) in the limit as long as (A1)-(A5) still hold. 



In order to avoid difficulties related to comparison results for sub and superso- 
lutions, we assume that F{x,u,p, M,l) = au for x large enough, say for |a;| > i?, 
where a is given by (A2). Taking into account the fact that the truncation on the 
distance function implies that d{x) = 1 for x large enough, the equation outside a 
large enough ball reduces to 



1 



+ -Y^{x) ■ Du^ - g] = 



which can be treated by a slight adaptation of the technics used in Section [2] as we 
will see it later on. For other extensions of F , additional conditions are typically 
needed to handle the growth (typically linear) of the solutions at infinity. 

Here it is unavoidable to impose additional assumptions on 7,5,^, to satisfy 
the integrability assumption (2.7), i.e. to balance the growth u(x + j(x, •)) with the 
decay of /i at infinity for solutions m of (1.1) and (5.4). We are going to use (BC3) 



and refer the reader to the discussion at the end of Section [2] 

We just recall that, in the case when g has compact support, the solutions are 



expected to be bounded by Lemma 2.1 and no additional assumption on j and /i 
is needed. On the contrary, if, for example, g = 1, then the integral of g in (2.1) 



suggests a linear growth and one has to impose suitable hypothesis on 7, j and /i in 
order to satisfy (2.7). Moreover, if we were considering more general extension of 



F, we would need a framework where we can compare sub and supersolutions with 
linear growth. Our restrictive extension allow us to avoid such (useless) technicali- 
ties. 

Theorem 5.4. Assume that (A1)-(A5) and (BC1)-(BC3) hold. Then, for any 
K > 0, there exists a unique continuous viscosity solution 0/ (5.4) which is 



uniformly locally hounded. Moreover, as k ^ 0, converges locally uniformly to 



the unique viscosity solution u of (1.1 1 -(1.2) 



In the proof we use the following lemma. 

Lemma 5.5. Assume (BC1)-(BC3). There exists a function 6 : 
that 

7(2;) • D9{x) > 1 for X in a neighborhood W of . 
Moreover 9 satisfies 

|6'(x)| < c(l + |a;|) inR^ , 

for some c > 0. 

We prove this result after the proof of Theorem |5.4[ 



such 
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Proof of Theorem \5.4\ We just sketch the proof of the existence and uniqueness 
of when g is not compactly supported. This case involves the function 9 of 
Lemma |5.5| while the other case is easier and involves a similarly defined but 
bounded function 6 (where D9 • 7 > 1 only on a compact set). 

The strong comparison principle (and hence uniqueness) for (5.4) holds by stan- 
dard argument and a slight modification of the argument of Section [2] that we 
explain now. If u is a subsolution of (5.4) then we have 



KaUn + 7(2:) • — g — 



in 



where R is defined above, Bp^ is the ball centered at with the (large) radius R. 
A slight modification of the arguments of Section j2j shows that, if y £ -B^ and if 
Xy{s) e i?| for s e [0,t) then 

Uniy) < / g{Xy{s))exp{-Kas)ds + Ui^{Xy{t))exp{-Kat). 

JQ 

Using this result, we can reduce to the case where the maximum points are in a 
fixed compact subsets of and then classical comparison arguments apply. 

Using Lemma 5.5 and (A2), it is easy to check that, choosing first C2 > 
and then Ci large enough, ±(Ci + C20{x)) are respectively viscosity super and 
subsolutions of (5.4). Then we can apply Perron's method to obtain the exisitence 
of a solution such that 

-(Ci + C20ix)) < u,{x) <Ci+ C20{x) in . 

Since the u^s are locally uniformly bounded, we can use the half-relaxed limits 
method. We rewrite (5.4) in the following equivalent way as G'„(x, m, Du, D^u,X[u]{x)) 
in where 



G^{x,r,p, X,l) 



F{x,r,p,XJ), 



F{x,r,p,X,l) +-f{x)p~ g, 



for X G Q,, 
for X € fj'^. 



G{x,r,p,XJ) = < 



As in the proof of Theorem |5.1[ we compute the half relaxed limits and find that 

F{x,r,p,X,l) when x E il, 

m.m{F{x,r,p, X,l),j{x)p — g} when x £ dft, 
7(x) ■ p — g when x € Cl'^^ 

and that G is like G with a max replacing the min, and we find that w is a viscosity 
subsolution of the G{x, u, Du, D'^u,2[u]) = and while w is a viscosity supersolution 
of the equation G{x,u, Du, D'^u,I[u]) = in i?^. We conclude as before that 
u — u—: u and ^ u locally uniformly. □ 

Now we give the proof of Lemma |5.5[ 



Proof of Lemma 5.5 This is a routine adaptation of classical arguments. Taking 
5 > small enough and denoting by Ds := {x e M^; d{x) < 6} where d is defined 
in Section [4T] we can solve the problem 

j{x) ■ Dw{x) — 2 \ii Ds , w — Q on dDg . 



(5.5) 



Indeed, arguing as in Lemma 2.1 with g ^ 2 , we have, for any y G Ds 

wiy) - 2t^ for = M{t > 0; Xy{t) e ODg}, 

and the function w is finite (thus well-defined) because of (BC2). 
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We prove that w is locally Lipschitz continous in Dg if 6 is so small that by 
(BCl), 

7(2;) • n{x) > — in As = {x : \d{x)\ < S}. 

We first check that w is Lipschitz continuous in A^. Let x,y ^ As, f{t) :— d{Xx{t + 
Ty)), and note that if > Ty, then 

fit) = X,{t + tD) ■ Dd{X,{t + r'y)) = j{X,{t + tD) ■ n{X,{t + r^)) 

for t e (0, Tx — Ty). We integrate from to — and use (BCl) to find that 



]\<\d{Xx{T'y))\<\Xx{T'y)-Xy{rl)\, 



(5.6) 



where the last inequality is a consequence of the definition of the distance of the 
point Xj.{Ty) to the boundary. Then if L is the Lipschitz constant of 7, inequality 
(4.11 1 holds and we may use e.g. (BC3) to obtain that 



<gLc(l + fl)|^. 



y\ 



(5.7) 



where R — max^^-^^ \x\. It follows that w is Lipschitz in As. 

Let X, y e Ds \ As be near one another and take a T > such that Xj.{T) E As. 
Such T exists and T < c(l + |a;|) by (BC3). By inequality (4.11 ), we can (and do) 
take y close enough to x so that also Xy{T) € As. Then — T 



Ty = T + ^rpy and hence by (BC3) and inequalities (5.7) and (4.11 1, 

^|rf - T^\ < e^^(i+^)|X.(T) - XyiT)\ < e^'^(i+^)e^^(i+l-l) |x - y\. 

This completes the proof of local Lipschitz continuity of w. 

The next step is to regularize w through a classical convolution argument to 
obtain the smooth function 0. But since w is only locally Lipschitz continuous, 
we have to regularize locally and use a covering argument to build the global reg- 
ularization of w. The covering argument is completely standard and will not be 
detailed here. 

Locally we define Wg{x) = w * Peix) for x G Ds where < e < | and p^{x) is 
the standard mollifier, i.e. a positive C°°-function with mass one and support in 

< £. By the regularity of w, Dw exists a.e. and hence equation (5.5) holds a.e. 
It follows that {Dw ■ ^) ^ — 2 m D s. By the definition of the convolution and of 
Pg, the Lipschitz continuity of 7, and the local boundedness of Dw, we are lead to 

Dw^ ■ "f[x) = [Dw ■ 7) * pe(a;) + j Dw{y) ■ {-f{y) - -/{x))ps{x - y) dy 
>2-\\Dw\\L^(^B{x,e))Lfe in D|. 
Hence for any bounded subset K C Ds we can take e — Ek so small that 

7 • Dw^^ > 1 in K. 

Finally, the bound on \9\ follows directly from a similar bound for w and a 
suitable (local) choice of e. The bound for zz; is a direct consequence of Assumption 
(BC3) and Lemma [2T| □ 
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